Anomalous Spin segregation in a weakly interacting two-component Fermi gas 
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We explain the spin segregation seen at Duke in a two-component gas of ^Li [3| as a mean-field 
effect describable via a coUisionless Boltzmann equation. As seen in experiments, we find that slight 
differences in the trapping potentials in the two spin states drive small spin currents. Hartree-Fock 
type interactions convert these currents into a redistribution of populations in energy space, and 
consequently a long lived spin texture develops. We explore the interaction strength dependence 
of these dynamics, finding nontrivial dependence on system parameters and close agreement with 
experiment. 



Spin waves in dilute cold gases such as spin polar- 
ized Hydrogen 0, Q or *'^Rb 0, i, d, 0] are driven 
by exchange effects and are manifestations of quantum 
coherence in a non-degenerate gas. Recent experiments 
at Duke [l[ on ^Li, with the interactions tuned to be 
very small, saw an unexpected spin segregation describ- 
able as a very long wavelength, low frequency, longitu- 
dinal spin wave. Their attempts to explain this behav- 
ior in terms of hydrodynamic spin-wave theory failed. 
Here we use a coUisionless Boltzmann equation to ex- 
plain this anomalous phenomenon. Our approach was 
motivated by discussions with John Thomas, who has 
recently explored a simplified version of this theory in a 
work with Du, Luo, and Clancy A concurrent study 
by Piechon, Fuchs, and Laloe [9(] reached similar conclu- 
sions. Spin waves in this coUisionless Knudsen regime 
have been studied extensively, both theoretically and ex- 
perimentally [HI, EE 111. 



In the Duke experiments []| a cold gas {T/Tp ~ 4) 
of roughly 2 x 10^ ^Li atoms, in their lowest hyper- 
fine state (denoted a were prepared in an optical 
plus magnetic trap with a trapping potential of the form 
Ui — ^mujjfr^ + ^muj^ ^Z'^, with = 27r x 4360Hz and 
ujz,i = 2ttx 145Hz. A radio pulse was used to coherently 
transfer atoms into a superposition of the | and the next 
hyperfine level (denoted t)- The subsequent dynamics 
were observed for several different bias magnetic fields, 
hence several different scattering lengths. 

When the scattering length was small and negative 
they observed that after ~ 100ms of evolution, the two 
components of the gas segregate axially with the | com- 
ponent moving inward, and the other moving outward. 
This spin texture persisted on timescales of a few seconds, 
much longer than the timescale for small oscillations. 
When the sign of the scattering length was changed, the 



"f moved outward and the j moved inward. When the 
scattering length was tuned to zero, no dynamics was 
seen. 

This behavior was surprising. According to hydrody- 
namic spin wave theory the characteristic timescale 
for any oscillations should be given by the axial oscilla- 
tor time, roughly two orders of magnitude faster than the 
observed dynamics. Moreover the only possible mecha- 
nism for driving this spin segregation is the very slight 
difference in the axial trap potential ^[(f/f ~ Ui)] ~ 
27r(4.4 X 10~"')Hz//im^. It is surprising that such a small 
difference in the trap leads to such dramatic density re- 
distributions. 

An intuitive picture of these dynamics is presented by 
Du et al. Q. They note that since the time-scale of 
spin rearrangement is long compared to the oscillation 
period, local equilibrium is never attained. Instead, each 
atom's spin dynamics is controlled by a mean field, av- 
eraged over its periodic trajectory. Low energy atoms 
that spend more time in high density regions experience 
a greater mean field. The [ atoms see a slightly weaker 
trapping potential and hence, for a given energy, have 
trajectories which extend over more space. This results 
in those atoms seeing a smaller mean field. The net re- 
sult of the subsequent dynamics is a spin segregation in 
energy space. 

Here we show how this behavior plays out in phase 
space. Following standard arguments [HI, [IBIj we derive 
an effective ID coUisionless Boltzmann equation. Solving 
this equation numerically, we reproduce the experimental 
observations. 

At the temperatures of interest (T ^ 27/iK) one only 
needs to consider s-wave scattering and the Hamiltonian 
reduces to 




where the field operators obey the fermionic equal time anti-commutation relations {'5^(ri, t), (r2, i)} = 
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FIG. 1: (Color Online) Time evolution of one-dimensional densities S| and S|: The up- (black, solid), and down- (red, dashed) 
spin density in space (horizontal axis in units of = 20 (moj^)"^^^) for (a) a — 4.5a_B, (b) a = and (c) a = — 4.5a_B, where 
as is the Bohr radius. Time runs from bottom to top, with each offset density profile is separated by 25ms. No dynamics are 
seen on this scale for a = 0. 



(5(ri — r2)(5^g.', and g = with s-wavc scattering 
length a. Near the magnetic fields of interest a{B) = 
-3.5(5 - Bo)aB/G We have set fi = 1 throughout, 
and we work in the Larmor frame rotating with a fre- 
quency equal to that of the transition for a uniform 
gas. 

Given the small scattering lengths, and low densities 
in this experiment (a ^ 4as, n lO^'^^cm"'^), the mean 
collision time r = l/(nav) ~ 10 s is much longer than 
the timescale of the experiment, and interactions only 
enter at the mean field level. From experimental studies 
of relaxation in a single component gas it appears 
that the time between background collisions Tb, due to an 
imperfect vacuum, is also on the order of several seconds. 
For times short compared to r and Tb, one can describe 
the system in terms of a collisionlcss Boltzmann equation, 
where the system is described by a Hartree-Fock approx- 
imation. The long timescales involved in coUisional re- 
laxation also explain why a simple hydrodynamic theory 
does not capture the physics of the phenomenon. 

Following [lit ^se the Heisenbcrg equations for 
'I'cr(r, t) to derive the equations of motion for the spin 
dependent Wigner function 

V = /TT(p,R,t) /u(p,R,t)\ 

f,,'{p,R,t) = y'dre-P-(vI/t(R-£,t)*^,(R+£,i)), 

which is the quantum analogue of the classical distri- 
bution function. Here p represents the momentum, 
r = ri — r2 is the relative coordinate and R = '"^'^'"^ is the 
center of mass coordinate. Treating interactions in the 
Hartree-Fock approximation [3| , the resulting equations 
involve densities and currents such as {s^ ^' (R, t)) = 
(vI/t(R,t)vI/^,(R,f)) = J^f^^,(^p'^n,t) and 
(j,,.(R,t))=/(^p/,,.(p,R,i). 
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FIG. 2: (Color Online) Top: Central spin density difference 
{sz{0,t) = S|(0,t) — S|(0,f)) normalized to the total central 
density (s = S|(0) -I- S|(0)) for various scattering lengths. 
The total central density is constant in time. Bottom to 
top: a = —4.55 as, a = —2.45 as, a — as, o, — 2.8 as, 
and a — 4.55 as- (b): First 100ms of evolution of s^/a for 
a = 4.55 as, and different values of 5u) = tj^ — lJ|: thin 
(red) - 27r X 0.15 mHz, green - 2n x 1.5 mHz, thick (blue) - 
27r X 15 mHz. The green curve corresponds to the experimen- 
tal value of Suj. 

We define s^a = Scr, s-^i = s+, s^f = s_ and analo- 
gously for the spin currents. Here s± refer the spin rais- 
ing and lowering operators that are related to the trans- 
verse components of the spin and Sy in the usual way 
.s± = Sx i iSy. Throughout this paper we use upper-case 
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ward scattering, while the off diagonal components rep- 
resent the exchange contribution. 

Given the differences between the radial and axial trap- 
ping frequencies, the relevant dynamics of the system 
are one-dimensional. In the nondegenerate limit that 
we consider here one can reduce (Eq. 3]) to a one di- 
mensional Boltzmann equation by making the ansatz: 
/(p,R,t) = f{p^,pz,R±,Z,t) = f{p^,R^)f{pz,Z,t), 
where the notation _L has been used to denote the trans- 
verse directions. The initial distribution is assumed to 
be a stationary state of the non-interacting Boltzmann 
equation for spin-down particles: 
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FIG. 3: (Color Online) Spin segregation in phase space. 
Lighter colors represent higher density. The total phase space 
distribution (left column), up-spin distribution (middle col- 
umn) and down spin distribution (right column) ai t — (top 
row), t = 100ms (center row) and t = 200ms (bottom row) 
for a = —4.5 as are shown. The phase space distribution is 
only a function of energy, but it is not of a simple Boltzmann 
form. 



letters to denote the components of position and momen- 
tum, and lower-case letters to denote the spin degrees of 
freedom. These transverse spin components represent a 
quantum coherence which is not captured by a classical 
model of a two component gas. The z-spin density is 
SzCRjt) ~ S|(R, t) — S|(R, t), and the particle number 
is TV = /rfR s(R,t) where s{R,t) = st(R,<) + Si(R,i). 
Assuming that all potentials are slowly varying in space 
and time, one finds 

O -I 

V+ P.Vr'f =z[V,'f] + -{VrV,Vp'f} (3) 

at ra I 



where the potential matrix is 



V = 



T 

gs- 



u 



(4) 



and the effective potentials arc U^'^^{R,t) = U-\{R) + 

gsi{Il,t) and Uf^'{R) = Ui{R) + gs|(R,t). Commu- 
tators and anti-commutators are respectively given by 
square brackets and braces. The diagonal terms of the 
potential matrix include the direct contribution to for- 



where the prefactor A = NlUzuJI / [ksTY is defined such 
that the constraint N = j (iR(s|(R, t) + sx(R, t)) is sat- 
isfied. 

By assuming the above Gaussian den- 
sity profile, we express the 3D density as 



„3D 



',(p,i?^,z,t) = 4/(|H^/, 



^ ^ ^^faa'iPz-,Z,t). Spatial averaging in the trans- 
verse direction renormalizes the coupling constant by 
i Finally, integrating the transverse momentum 

coordinates yields an effective interaction strength 
9^° = ^gN{muj,){LOrlkBTf. 

Results and Discussion:— Working in units where 
lengths and momenta are measured in terms of the os- 
cillator length {mu}z)~'^/'^ and {mkBTy^^ respectively, 
we use a phase space conserving split-step approach |17| 
to integrate Eq. (jH), using a 20 by 20 by 800 grid in 
Z-pz-t with step sizes 6Z = 20 (mwz)"^''^, Spz = 
0.6 [mkBTY''^, and 5t = 0.04^. We verified that our 
grid was sufficiently fine so that our results no longer 
depended on the step sizes. 

In Fig. [1] we show the time evolution of the density 
profile for the up (black/solid curve) and down spins 
(red/dashed curve) in space for the first 200 ms, initial- 
izing the cloud in an Sx state. We plot the behavior 
for three different values of the scattering length, finding 
that both the magnitude and timescales are in agreement 
with experiment [l|. 

In Fig. [2l^a) , we quantify the magnitude of the spin seg- 
regation by plotting the central density difference (s^ = 
S| — S|), normalized to the total density (s = st + ■^i) '"^ 
a function of time for a range of scattering lengths. This 
quantity peaks near t ^ 200 ms. We extract the timescale 
associated with spin segregation by taking the slope of 
the graph at small times (Fig. Wih) green curve) to find 
_ 1/200 (ms)-i for a = 4.55aB . Furthermore, 
the figure reveals oscillations in the spin density differ- 
ence at a frequency ~ 27r300 Hz (~ ^lUz) that is weakly 
dependent on the interaction strength, corresponding to 
the lowest breathing mode of a two component Fermi 
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gas [18|. The amplitude of these oscillations depends on 
the difference in the trap frequencies seen by the t and | 
atoms. These oscillations are not captured in the analysis 
presented in Q. 

In previous experiments on bosons 0] , the spin dynam- 
ics were much faster than such collective modes. This dif- 
ference can be attributed to the ratios of the mean field 
interaction energy to the trap frequency A = g^^ /ujz- In 
the current experiment A ~ 0.2 while in A ~ 10. 

Fig-Etb) shows that both the magnitude and timescale 
for spin segregation seen in [l| is strongly dependent on 
the difference in trapping frequencies {5ui ~ ~ ldi)- 
Had this frequency difference in [H been an order of mag- 
nitude larger (blue/ thick curve in Fig. [Jib)), the dynam- 
ics would have been much more complicated and much 
less dramatic. 

As previously discussed, an important observation in 
[1] was that the spin segregation in [l[ can be viewed as 
a segregation in energy space. We illustrate this effect 
in Fig. [3] by plotting the phase space distributions for 
a = — 4.5a_B for t ~ 0, 100 and 200 ms respectively. One 
sees that the phase space distributions are not separately 
functions of Z and pz, but instead depend on uj^Z"^ -f- 
p|/m. 

Finally we note that this spin segregation is very ro- 
bust. We can illustrate this by exciting a large amplitude 
spin dipole mode at i = 0. We find that spin segrega- 
tion occurs even as the t and | atoms slosh around in 
the trap, out of phase with one another. As may be 
expected for a gas in the Knudsen regime, oscillations 
on timescales much shorter than the interactions do not 
change the long term dynamics. Nonetheless, it would 
be interesting to observe this stability experimentally. 

Summary and Conclusions:— Using standard tech- 
niques [l6l |. we have derived a coUisionless Boltzmann 
equation which reproduces the anomalous spin waves 
seen in [l|. This is an exciting regime for spin waves, 
as the system is far from local equilibrium. Remarkably 
we find an ergodicity where the phase space distribution 
function is only a function of energy - but is not a simple 
exponential. 

Our numerical simulations indicate that this spin seg- 
regation depends strongly on the difference in the trap- 
ping frequencies seen by the two spin species. More- 



over, despite being in a nondegenerate regime, substan- 
tial quantum coherences are found in this system. We be- 
lieve that much can be learned from studying how these 
coUisionless dynamics evolve into hydrodynamics as the 
scattering length is made larger. 
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